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Abstract 


In this thesis, a robust optimal controller for the launch vehicle is designed. Also, per- 
formance of a hybrid closed loop system containing continuous time plant and discrete 
time controller with and without lifting is considered. Two cases are considered: 

case I : Controller designed for discrete-time model of the plant. 

case II : Controller obtained by descretizing the continuous time controller. 

The controllers are designed using the reduced order model of the launch vehicle. 
Controller order is further reduced by balanced realization truncation and frequency 
weighted balanced realization truncation techniques. A comparison is made of the per- 
formance of the system for different controllers with and without lifting the system. Also, 
a comparison is made of the performance of the various reduced order controllers with 
the previously designed controller for the same system. The comparison clearly in- 
dicates the benefits of lifting the system. It is found that performance with the 2 nd 
order controller in frequency weighted balanced realization for case II is better than 
any other controller in two different realization presented for different order. Bounds 
on allowed perturbations in system matrices are determined. Bounds are calculated 
with and without lifting the system. Larger perturbation bounds are obtained with 
lifting than without lifting. 


iii 



Dedicated To 
My Family Members 



Acknowledgements 


I take this opportunity to express my sincere gratitude towards my thesis supervisor, 
Dr. K. E. Hole for his constant guidance and encouragement throughout the prepara- 
tion of this thesis. Without his help, this work would never have been completed. 

I would like to thank all my friends for their good company. I am thankful to 
my husband, Rajeev, without whose loving attention and support. I would not have 
achieved this success. 

Lastly, words do not suffice to express my indebtedness and the deep affection I 
have for my parents and all my family members, for they have been my source of 
inspiration and moral' support. 


v 



Contents 


1 INTRODUCTION 1 

1.1 General 1 

1.2 Review Of The Existing Work 2 

1.3 Proposed Work In Brief 6 

1.4 Outline Of The Thesis 6 

2 FAST SAMPLING AND LIFTING TECHNIQUE 8 

2.1 Introduction 8 

2.2 Lifting Technique 9 

2.3 State Models And Transfer Functions 11 

2.4 Conclusions 14 

3 CONTROLLER REDUCTION 15 

3.1 Introduction 15 

3.2 Balanced Realization Truncation 16 

J 

3.2.1 Procedure For Balanced Realization Truncation 17 

3.3 Frequency Weighted Balanced Realization 

Technique 19 

3.3.1 A State Space Representation For W(z) and V(z ) 21 

3.4 Construction Of Balancing Transformation 23 

3.5 Implementation of the Controller Order Reduction Procedure 24 

3.6 Conclusions ' 27 


vi 



vu 

4 NUMERICAL EXAMPLE 28 

4.1 Introduction 28 

4.2 A Sufficient Condition For Stability Robustness 28 

4.3 Design Of Autopilot For Launch Vehicle 30 

4.3.1 System Model 30 

4.3.2 Robust Optimal Controller Design 31 

4.4 Conclusions 38 

5 CONCLUSIONS 46 

APPENDIX A 48 

APPENDIX B 63 


References 


64 



List of Figures 


2.1 The Closed Loop System 9 

2.2 The Fast Sampled Closed Loop System 10 

2.3 Blocking or Lifting on Serial to parallel conversion 11 

2.4 The Lifted Closed Loop System 14 

3.1 Redrawing of closed loop system with controller K r (z) 16 

3.2 Block diagram for obtaining a state space representation for W(z) ... 21 

3.3 Block diagram for obtaining a state space representation for V(z) ... 23 

4.1 Step-Responses of the 9 th order LQGd , 4 tA order LQGbd and 7 th order 

LQGfbd controllers on full-order (23rd) nominal system 39 

4.2 Step-Responses of the 9 th order LQGs, 3 rd order LQGbs and 2 nd order 

. LQGfbs controllers on full-order(23rd) nominal system 40 

4.3 Step-Responses of the 9 th order LQGd, 4 th order LQGbd and 7 th order 

LQGfbd controllers on full-order(23rd) perturbed system set-1 41 

4.4 Step- Responses of the 9 th order LQGs , 3 rd order LQGbs and 2 nd order 

LQGfbs controllers on full-order (23rd) perturbed system set-1 42 

4.5 Step-Responses of the 9 th order LQGd , 4 t/l order LQGbd and 7 th order 

LQGfbd controllers on full-order(23rd) perturbed system set-II 43 

4.6 Step-Responses of the 9 th order LQG S , 3 rd order LQGbs and 2 nd order 

LQGfbs controllers on full-order(23rd) perturbed system set-II 44 

4.7 Bode magnitude plot with 2 nd order LQGfbs controller 45 


vm 



List of Tables 


4.1 % overshoot with lifting 34 

4.2 % overshoot without lifting 35 

4.3 Bandwidth, Gain margin and phase margin with different controllers . 36 

4.4 Perturbation bound with different controllers 38' 

A.l Zeros of the continuous time full order system model 48 

A. 2 Zeros of the discrete time full order system model (sampling period 

=0.02 secs.) 4.9 


rx 



Chapter 1 


INTRODUCTION 


1 . 1 General 

Due to rapid development in the technology of digital computers and microprocessors, 
considerable attention has been focused on the study of digital and sampled-data con- 
trol systems. For such systems to obtain the best possible performance, it is often 
necessary to sample the signals for the various sensors and control effectors at different 
rates. A sampled-data system is considered a multirate system if sampling at different 
locations occurs at different rate. 

The two primary operations in multirate processing are decimation and interpola- 
tion. They enable the data rate to be altered without significant, undesirable effects 
of errors such as quantization and aliasing. Decimation reduces the sampling rate ef- 
fectively. Interpolation on the other hand increses the sampling rate. Fast sampling 
of the system results in periodic, discrete-time, linear and time-varying system. Fast 
sampling of the system at a multiple of the sampling frequency followed by lifting al- 
lows capturing of the systems intersample behaviour and yields a time invariant single 
rate system. This then permits standard order reduction ideas to be applied. Lifting 
technique is given in [1]. 

However, in [2], [3] general multirate sampled-data system is considered. Here 
all the inputs and outputs sampled at different rates are considered. It is shown that 
multiinput, multioutput, multirate sampled-data systems, while not generally periodic, 
are shift varying and, for any given multirate system, obtains an explicit expression for 
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a state space realization of its shift invariant equivalent. 

Multirate sampling schemes becomes necessary for systems with special data trans- 
mission links or special sensors and actuators. They are useful for improving the 
response of systems in which measurements are obtained at slow rate e.g. when lab in- 
struments are used. Multirate systems allows better control of what happens between 
sampling instants. In multivariable systems it may also be advantageous to have differ- 
ent sampling rates in different loops to reduce the computational load and to improve 
the numeric conditioning. Use of multirate sampling is also natural in multiprocessor 
systems [4]. High quality data acquisition and storage system takes advantage of multi- 
rate technique to avoid the use of expensive anti-aliasing analogue filters and to handle 
efficiently signals of different bandwidths which require different sampling frequencies 
[5], 

Simple linear controllers are normally preferred over complex linear controllers. For 
this reason, there is a desire to have methods available for designing low order con- 
trollers for high order plants. Such methods can be divided into two classes, direct, in 
which the parameters defining a low order controllers are computed by some optimiza- 
tion or other procedure; and indirect, in which a high order controller is first found, 
and then a procedure used to simplify it [6]. 

Any mathematical model, no matter how accurate, can only approximate the true 
behaviour of the system. This is because of various reasons such as linearization of 
nonlinear system about an operating point, order reduction and neglected coupling 
terms. Further, the parameters involved in such a description are often subjected 
to variations and uncertainties for several reasons such as environmental condition, 
ageing, wearing and manufacturing tolerance. It is therefore important to assess the 
robustness of the system stability to parameter variation. 

1.2 Review Of The Existing Work 

Although study of systems containing multiple samplers or hold operating at different 
rates began in the 1950’s, it has received renewed attention. Interest is partly due to 
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the importance of such control systems in practice since a very large number of ” real- 
world” digital control systems are multirate. Sklansky and Ragazzini [7], developed 
the frequency decomposition method for determining the transfer function of multirate 
system. The switch decomposition method, developed by Kranc [8], is an extension 
of the frequency decomposition method which allows rational sampling rate ratios. 
Ragazzini and Franklin [9], includes an excellent review of the frequency decomposition 
and switch decomposition methods. 

Kalman and Bertram [10], presented a method for forming a discrete time state 
model of an multirate (MR) sampled data plant. A frequency domain approach for 
discrete-time periodic time- varying systems appears in the work of Davis [11], Jury and 
Mullin [12], and Meyer and Burrus [13]. Nagy and Foias [14], have shown ’categorical’ 
equivalence between periodic discrete-time systems and certain kinds of linear time 
invariant (LTI) systems. 

Based on the results of [14] and [15], Khargonekar et al. [1], derived transfer function 
for periodic, discrete-time, linear, time- varying systems. It is shown that lifting has the 
effect of removing blocking zeros. They, then, applied this observation to several robust 
control problems, proving in particular that, for a single input single output (SISO) 
bicausal plant, a periodic controller can be designed so that the resulting system has 
an arbitrarily large gain margin. In [16], [17] it is shown that there is a limit to the 
gain margin achievable by a time invariant controller if the plant has an unstable pole 
and zero. 

To analyze the behaviour of continuous time periodic system Bamieh et al. [18], used 
a lifting technique similar to that used for discrete time periodic systems in [1]. They 
have shown there that lifting technique is applicable to all norm based optimization 
problems, and in particular to sampled-data versions of the quadratic regulator and 
optimal filtering problem. Yamamoto [19], also used lifting for sampled-data systems, 
but he lifted the state as well as the input and output. 

Many researchers have considered the control system design problem for multi rate 
sampled data plant (MRSDP) under various assumptions on the sampling rates and 
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the type of control used. In [20] authors analyze multirate pole placement assuming 
that the inputs and outputs can be paired off into loops having the same rates. Suc- 
cessive loop closure [21] is a heuristic technique applicable under the same assumption. 
An Linear Quadratic Gaussian (LQG) design method for pairwise single rate plant is 
reported in [22] and uses a mechanism from [10]. 

In [23], Franklin and Rahmani introduced a new optimal multirate control of linear 
continuous-time periodic systems. The control law is based on the continuous time 
linear quadratic regulator (LQR) problem. With the input being constrained to a 
certain piecewise constant signal. It is shown that multirate control scheme, when 
applied to LTI systems, can produce better response characteristics and a less LQR 
cost than that of single rate sampling if the state is sampled at the same rate in both 
the schemes. 

Another simplifying assumption is that all the inputs (respectively, output) are 
updated at a single rate. Another is so called ’’dual rate sampling” where all inputs 
are updated at one common rate and all outputs are updated at another common 
rate. State feedback problems under these frameworks are addressed in [24] and pole 
placement type design methods are given in [25] and [26]. 

In [3], Meyer studied the multirate LQG problem using standard LQG theory and 
without simplifying assumptions on the sampling. It is shown that multirate LQG 
problem can be solved without using periodic system theory or solving periodic riccati 
equations. Here multirate sampled data LQG problem is translated into an equivalent, 
modified, single-rate, shift-invariant problem via a lifting isomorphism approach. In [2], 
A new class of linear shift- varying operators that generalizes the notion of N-periodicity 
is defined. The stability theory for linear time-invariant plants with periodic digital 
controllers is considered in [27]. Here it is shown that gain margin can be arbitrarily 
assigned by suitable choice of sampling time and digital controller. 

Examples of direct methods for controller reduction include the work of Gangsaas 
et al. [28], Hyland and Bernstein [29]. For indirect method good number of references 
can be found in [6]. A sampled data controller reduction procedure is considered in 
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[30]. Here weighted balanced truncation is used to reduce the order of the controller 
designed for multirate sampled system using lifting technique. Anderson and colleague 
have developed a frequency weighted hankel norm technique for controller reduction 
[6], [31]. A frequency error bound for the frequency weighted controller order reduction 
when stability is preserved can be found in [32]. 

There are now at least three important and also rather popular state space based 
model reduction techniques, namely, truncation of internally balanced realization [33], 
[34] Hankel norm optimal approximation [35], [36] and q-covariance equivalent real- 
ization (q-cover) [37], [38], [39]. Other approaches employing balancing to achieve 
controller reduction includes work of Yousuff and Skelton [40] and Davis and Skelton 
[41]. A general information on controller design is given in [42], In [43] optimal hankel 
model reduction for nonminimal systems is discussed. Balancing transformations for 
unstable nonminimal linear system is given in [44]. 

A discrete time controller is obtained by finding the zero order hold equivalent of the 
open loop controller. It is shown in [45] that this may lead to poor performance because 
it does not seek directly to have the continuous time closed loop closly approximated by 
the sampled data closed loop. Stability of sampled data feedback systems is discussed 
in [46]. 

The stability robustness analysis in face of structured uncertainty for continuous 
time case had been given by Yedavalli [47], Barmish [48], Zhou and Khargonekar [49], 
Hole [50], Kosmidou [51]. A simple test for determining the stability of a characteristic 
polynomial whose coefficients vary around some set of nominal values has also been 
developed. For such type of interval polynomial, Kharitonov showed that stability can 
be determined from the stability of four fixed polynomials. Extension of Kharitonov 
theorem to discrete case are given in [52], [53]. 

In [54], smallest destabilizing perturbations for linear systems(both for continuous 
and discrete time case) have been derived. A sufficient condition for the stability 
robustness of linear systems with time varying structured uncertainty is given by Sobel 
et al. [55]. Here it is shown that the nominal eigenvalues lie to the left of a vertical 
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line in the complex plane which is determined by a norm involving the structure of 
the uncertainty and nominal closed loop eigenvector matrix. In [56], element bound 
for allowed perturbations in system matrices is derived based on discrete lyapunov 
stability approach. Niu et al. [57], have given bounds for allowed perturbations in 
system matrices using discrete lyapunov stability approach. 

Luo and Johnson [58], have given a saddle point inequality of the closed loop system 
performance to cope with noise uncertainty, and a sufficient condition for system sta- 
bility robustness under plant perturbation. These results are extended in [59], where it 
is shown that the new robust stability sufficient condition may become necessary and 
sufficient condition of system stability in certain cases. 

1.3 Proposed Work In Brief 

In the present work closed loop performance of a sampled data system via lifting tech- 
mque is considered. Based on results given in [1] lifting technique is used to replace the 
periodically time-varying sampled-data system by a single-rate time-invariant system. 

Based on result given in [6] controller order reduction by balanced realization trun- 
cation and frequency weighted balanced realization truncation (both for discrete time 
controller and controller obtained after discretizing the continuous time controller) is 
considered. Also bounds on maximum perturbations in entries of system matrices are 
calculated using results given in [59]. 

The effectiveness of the theory has been shown by designing digital autopilot for 
launch vehicle. 

1.4 Outline Of The Thesis 

The Chapter 2, presents fast sampling followed by lifting for the multirate hybrid 
closed-loop system containing plant and controller. A time-invariant system results. In 
Chapter 3, Controller order reduction by balanced realization truncation method and 
frequency weighted balanced truncation method is considered. In case of frequency 
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weighted balanced truncation method, weighting functions for preserving the closed 
loop transfer function is presented. A procedure for controller order reduction in both 
the cases is given. In Chapter 4, stability robustness of the systems with structured 
time varying uncertainties is considered. Also a practical example of designing digital 
autopilot for launch vehicle is given. The Chapter 5, concludes the thesis. 



Chapter 2 


FAST SAMPLING AND LIFTING TECHNIQUE 

2.1 Introduction 

In this chapter, the fast sampling and lifting technique is discussed. The periodically 
time-varying sampled-data system is replaced by a single-rate time-invariant system. 

The application of the ordinary z-transform to the analysis and synthesis of sampled 
data system gives information which can be used to evaluate the variables at sampling 
instants only. However, the behaviour of the system between sampling instants is 
of considerable importance. For instance, system may have considerable overshoots 
between sampling instants during the transient period, even though the response at 
sampling instants has settled completely. 

To take into account intersampler behaviour of the system, multirate sampling 
as described in [30] is used. Then the system becomes periodically, linear, time-varying 
sampled-data system. Using lifting technique as described in [1], this system is replaced 
by a single-rate time-invariant system. 

Fast sampling and lifting technique is described in Section 2.2. Also, the state 
model of the lifted plant and lifted controller are presented. 


8 



9 


2.2 Lifting Technique 

In this section, the problem of replacing the periodically time-varying sampled-data 
system by a discrete-time, linear, time-invariant system is considered. 

Consider a hybrid closed loop system where the plant is continuous-time and the 
controller is discrete-time (controller may be designed for the discrete time model of 
the plant, or, may be obtained by descretizing the continuous-time controller). 

This closed loop system is shown in Figure 2.1. 


Sampler Controller Hold Plant 



Figure 2.1: The Closed Loop System 


P stands for the p x m continuous-time plant, 

K for the m x p discrete-time controller, 

S for the sampler with the sampling period r, and 
H for the hold element, here assumed to be a zero order hold, 
m for input of the plant, 
p for output of the plant. 

This is a periodically time-varying sampled-data system.The aim is to replace the 
system in Figure 2.1 by a time-invariant system. 

To take into account intersample behaviour of the system in Figure 2.1 fast sampling 
of the system with the sampling time r/N is introduced. Here r is the controller 
sampling time and N is a positive integer, and r/N is chosen to be a submultiple N of 
r. 
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The fast sampled closed loop system is shown in Figure 2.2 



Fast Fast 

Hold Sampler 


Figure 2.2: The Fast Sampled Closed Loop System 

The sampling rate can be chosen based on criteria given in [4] [60]. Normally, t/N 
is chosen to be much smaller than the smallest time constant in the scheme of Figure 
2.1. So the sampling rate can be chosen to be smaller than the inverse of 5 x closed 
loop bandwidth. 

The sampled system in Figure 2.2 is a multirate N-periodic discrete-time system. 
Lifting involves passing from an N-periodic linear p x m discrete-time system to an 
equivalent pN x mN discrete-time linear time- invariant system. Lifting is also known 
as ’’Blocking” or ’’serial to parallel conversion” [42], Using this technique the multirate 
system is changed into a single rate system. This technique can be easily understood 
by looking at Figure 2.3. 

In Figure 2.3, the upper system has signals arriving say at 1/secs., while the lower 
system has one 3 vector arriving every 3 secs.. Likewise, the upper system has one 
output coming out every second, and the lower system has one 3 vector arriving every 
3 secs.. The total information entering and leaving the two systems is identical; it is 
just presented differently. This shows that a fast system is slowed down and is replaced 
by an equivalent slow system [42]. If system matrices of a state- variable description 
of the fast system are known, the corresponding matrices for the slow system can be 
easily found [1]. This is presented in next section. 

Here the equivalence between periodic discrete time system and linear time invariant 
system is an isomorphism of the systems in the sense that both essential algebraic and 
analytic properties of the systems are preserved. In particular, the lifted system is 
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Figure 2.3: Blocking or Lifting on Serial to parallel conversion 


stable if and only if the N-periodic system is stable, and in this case the operator norms 
(associated with regarding the system as an operator mapping square summable input 
to square summable output) are equal. These results are given in [1], [27]. 

To obtain a time-invariant system, we lift the system in Figure 2.2. 

Transfer function of the lifted system is derived next. 

2.3 State Models And Transfer Functions 

Consider the following continuous-time model of a linear time-invariant plant P(s) : 

= A p x r (t) + B p u(t) 
at 

z p (t) = Mp Xp (■ t ) + Vp u {t) 

y p (t) = C p Xp(t) + D p u(t) (2.1) 

where x p is the state vector of dimension (n p x 1), u is the control vector of dimension 
(n u x 1), y p is the output vector of dimensions (n y x 1), and z p is the measurement 
vector of dimensions (n z x 1). 
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The discrete-time approximation of fast sampled system is obtained from (2.1) for 
ampling period of ( t/N ) secs., and a zero order hold. Where r is controller sampling 
time and N is a positive integer. 

The discrete time model thus obtained is described by the following equation: 

x d (k + 1) = A d x d (k) + B d u(k) 
z d (k + 1) = M d x d {k) + V d u{k) 
y d (k) = C d x d {k) + D d u(k) 


where (2.2) 

A d = exp(A p r/N) 


B d = 

rr/N 

j exp(A p t)dt Bp 

Jo 


c d = 

c p 


D d - 

D p . 

(2.3) 


The dynamic controller K c (t) for the continuous-time model of the plant is described 
by tTe following equation: 

ffegM = A c x c (t) + B c z p (t) 

dt 

u(t ) = C c x c (t) + D c z p (t). (2.4) 

whore x f . is the controller state vector of dimension (n c xl ),u and z p are the control 
and measurement vectors, described in the plant model. 

The discrete-time controller can be found by descretizing the continuous- time con- 
troller given in equation (2.4), or, can be found for the discrete-time model of the plant 
given in equation (2.1). In both the cases the sampling time is r secs. 

For convenience following notations are used. 

LQG s l controller obtained after descretizing the continuous-time controller, 

LQG)£ controller obtained for the discrete-time model of the plant. 

The discreie- time controller K(z) is described by the following equation: 


x(k + l) = Ax(k) + B z d (k) 
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u(k) = C x{k) + D z d (k) 


(2.5) 


where 3? is the controller state vector of dimensions (n x 1) , u and z d are the control 
and measurement vectors, described in the plant model(2.2). 

Now the plant described by equation (2.2) is N-periodic discrete-time system. To 
replace this system by single-rate, linear, time-invariant system, the system is lifted as 
described in [lj. Then the state model of the lifted plant p is 


X L (k + 1) = A L X L (k) + B L U(k) 
Y L (k + 1) = C L X L (k) + D L U(k) 


where 


Al = A N d 
B l = [A N d ~ l B d . 
C L = [CjA T d Cj 

D d 

C d B d 


A d B d B d ) 

■ (A T d ) N ~ l C T d ] T 
0 

D d 


D l = 

C d A d ~ 2 B d C d A d ~ 3 B d 

The realization (2.7) is minimal if (2.2) is minimal. 

The lifted controller « is given as below 

k(z)=E 1 K(z)E 2 

where 

E\ — ( I'm Im ■ • • Ira ) ^ 

E2 = ( Ip Op - . • Op ) £ ffi-pxpN 


0 

0 

D d 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


(2.9) 


I n l n x n identity matrix. 

Op l p x p zero matrix. 

E 2 in equation (2.9) corresponds to a slow (every r seconds) sampler, which passes 
through only the first element of an input vector anH is in the off mode when the 
following N-l elements of the input vector arrive. 
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E\ corresponds to a r-second zero order hold. 

Ei and E% in the single-rate system correspond to interpolator and decimator in the 
multi-rate system. Introducing 

P = pEi ( 2 . 10 ) 

the periodically time varying sampled-data system in Figure 2.1 can be replaced by 
a time invariant system in Figure 2.4. 


Controller Plant 



Figure 2.4: The Lifted Closed Loop System 


2.4 Conclusions 

In this chapter, the lifting technique to replace the periodically time- varying sampled- 
data system by a single-rate, discrete-time, linear, time-invariant system has been 
presented. Also, state models of the lifted systems are presented. 






Chapter 3 


CONTROLLER REDUCTION 


3.1 Introduction 

The controller design methods specially LQG (Linear quadratic gaussian) or H„ for 
physical systems with high order models normally result in a controller of order equal to 
that of model. Low-order controllers are normally preferred to high-order controllers, 
given comparable performance: there are fewer things to go wrong in the hardware; 
their operation is easier to grasp at the conceptual level; that is, one is more likely to be 
able to identify parts of the controller as achieving certain subgoals of control; and m a 
discrete-time implementation the computational requirements are less. Therefore, it is 
desirable to reduce the controller order. Here indirect approach to low order controller 
design is considered. This uses LQG method to design a high order controller and then 

this controller is approximated by a low order contr 

In this chapter, the procedure for carrying out this approximation is considered. In 
Section 3.2, the balanced realisation truncation technique is presented and in Sect.on 
3.3, the frequency weighted balanced truncation technique for controller order reduction 
is detailed along with a weighting function for preserving the closed loop transfer 
function. Also, a balancing transformation matrix is presented. 
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3.2 Balanced Realization Truncation 

let Piz) and K(z) be the tranfer function matrices of lifted plant and controller as 
described in equations (2.10) and (2.5) respectively. 

Here Kiz) is a stabilizing high order dynamic controller with unity negative feed- 
back. Let K r (z) be a low order controller, which we are seeking. Closed loop system 
with controller K r (z) replacing X (z) is shown in Figure 3.1. 



Figure 3.1: Redrawing of closed loop system with controller I\ r {z) 


Then we have the following problem: 

Given transfer function matrix K(z) of order n with all poles located inside the unit 
circle in the z-plane centered at origin, find K r {z) of order r < n with all poles located 
inside the unit circle in the z-plane centered at origin, and in addition such that 

J = II K{z)-Kr(z) Hoc (3.1) 


is minimum. 

The notation 11.4(2)1100 means sup cWr [-4(a)]- Here, a max [-4(2)] denotes the maximum 

■0[0,2n] 

singular value of matrix -4(z) and z — e je . 

Condition (3.1) and the condition that K and K r have the same number of unstable 
poles suggest the following procedure for constructing reduced order controller.Let 


K{z) = K + (z) + K-{z). 


(3.2) 
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where K+{z) is strictly proper, with all poles outside the unit circle and K-(z ) has all 
poles inside the unit circle. Now the reduced order controller K r is represented as 

K r (z) = K + (z) + K- r (z). (3-3) 

where unstable part of K(z) is copied into K r (z) and K- r (z) has all poles inside the 
unit circle. And A_ r (^) is chosen to minimize J over all K- r {z ) of prescribed degree. 
Here the requirements of K(z ) and K r (z) to have same number of unstable poles and 
minimization of criterion given in equation 3.1 are sufficient conditions, not necessary 
conditions, for stability with the reduced order controllers [60], 

To solve this problem, at least in an approximate way, the balanced realization 
truncation technique [6], [60] is applied to stable part of K. 

3.2.1 Procedure For Balanced Realization Truncation 

Given an n-th order, linear, time invariant, and asymptotically stable controller with 
transfer function matrix K(z), where 

K(z ) = C{zl-A)~ l B (3.4) 

with {/l, B, (7} minimal. Let U, Y be the solution of the following Lyapunov equations 

AUA t + BB t = U 

A t YA + C t C = Y (3.5) 

when 1 U and Y are the infinite time controlability and observability Gramians. Con- 
sider a coordinate transformation matrix T. Let 

A b = T~ l AT 

B b = T~ l B 

C b = CT. (3.6) 

such that 

A 6 S Af + B b B b = S 
A^ S Aft + C b C), = E. 


(3.7) 
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where 


E = diag(a 1 ,a 2 , ... ,a n ) a x > a l+i 


Then the matrix T is called a balancing transformation and {At, Bt , Ct } is a balanced 
realization of {A, B, C}. The matrix E is both the controllability and observability 
Gramians. The Hankel singular value a t are obtained as positive square roots of the 
eigenvalue of UY . 

Now, for controller reduction partition the system {At, Bt, Cb} as 


At = 


A n 
A \ 21 



B h 




(3.8) 


and set 

A r — An ) B r = B x , C r — C\. (3-9) 


Then the reduced order system {A r ,B r ,C r } is a good approximation of system 
{A,B,C} ifcr r >or r+1 . 

In fact we have the following two properties 

• Subsystems {A li ,B i ,Ci}, i=l,2 are asymptotically stable if ay > oy + i [34]. 

• There exists a frequency domain error bound for the balancing approximation 
[36] 

|| C ( zl — A ) -1 B — C r (zl — A r ) 1 B r ||oo< 

2 ( oy+i + + <y n ) (3.10) 


Now if cx t+2 , . . . , cr n are much smaller than oy+i, from (3.10) it is clear that bal- 
anced realization truncation necessarily comes close to the optimum. Comparisons on 
practical examples suggests that balanced realization has its own in built frequency 
shaping which appears to enhance more often than degrade a controller approxima- 
tion. This reason and its relative simplicity are behind our preferance for the balanced 

realization approach. 
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3.3 Frequency Weighted Balanced Realization 
Technique 

In this Section the order of the controller K of the setup in Figure 2.4 using the approach 
given in [30] is considered. 

Any controller reduction procedure should take into account the existence of the 
plant. Controller reduction should preserve closed-loop stability, and as far as possible, 
the closed loop performance and closed loop transfer function. This generates frequency 
weighted approximation problems. The choice of frequency weight is influenced by the 
choice of criterion used [6]. Here the aim of controller reduction is the preservation of 
the closed loop transfer function. 

The closed loop transfer function matrices with K(z), K r (z)( where K{z), K r (z) 
are as described in Section 3.2 ) are 

T(z) = P(z)K(z)E 2 [I + P(z)K(z)E 2 }~ 1 

%(z) = P(z)K r (z)E 2 [I + P(z)K T (z)E 2 }- 1 . (3.11) 

Where P(z), E 2 are as described in equation (2.9-2.10). Consider the difference 

T(z)-T r (z) = 

P(z) K(z) E 2 [I + P{z) K(z ) E 2 ]- x - P(z) K r (z ) E 2 [ I + P(z) K r (z ) E 2 j" 1 . (3.12) 

To a first order approximation in K(z) — K r (z). there holds 
T(z) - T r (z) » 

[ I + P(z) K(z) E 2 3- 1 P{z) [ K(z) - K r {z) ]E 2 [I + P(z) K(z) E 2 J" 1 . (3.13) 

And this suggests the following approximation problem, Find I\ r of nominated degree 
such that 

K(z) and K r (z) have same number of unstable poles, and 

]| W{z) [K(z) - K r (z)]V(z) Hoc is minimized. 

IF(-) and V(z) obtained from eqn. (3.13) are given below: 

W(z) = [I + P(z)K(z)E 2 }- 1 P 

V(z) = E 2 [I + P(z)K(z)E 2 ]-\ 


(3.14) 
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are the weighting functions. 

o solve this problem, at least in an approximate way, the frequency weighted 
balanced truncation technique [30] is applied to the stable part of K{z) (As mentioned 
in Section 3.2, the unstable part of K(z) is copied into K r (z ) ). 

Consider asymptotically stable frequency weighting functions (3.14) and associated 
minimal state variable realizations: 


W ( z ) = C w [zl -A w ]~ l B w +D w 

V(z) = C v [zl - A v ]~ l B v + D v . (3.15) 


II (~) and V (z) arc stable when the closed loop 7~(z ) is stable. 

IIoic the basic idea is to change the controllability and observability gramians to 
teflect the intioduction of the frequency weighting. The frequency weighted transfer 
function II (, 2 ) Ii(z) V (z) has a representation with the following state space matrices 


Let 


and 



» 

B W C B W DC V 


B W D D v 

A = 

0 

A 

BCy 

B 

~ 

BD V 


0 

0 

Ay 


B v 

ry 

c w 

D W C 

D w D Cy 

m 






u n Un 

U n 




u = 

v& U 

co 





B\z U23 

u v 



(3-16) 


(3.17) 


Y = 


Y w Y a Y r , 
Y? 2 Y V 2! 

\ r T yT v" 
M3 M3 x v 


be the solution of the following lyapunov equations 


AUA t + BB t = u. 
A t YA+C t Q = 


(3.18) 


(3.19) 
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Now U and Y can be regarded as the frequency weighted controllability and ob- 
servability gramians for the original controller K(z). 

Consider a coordinate basis change to {A, B , C , } which makes 

Unew = Ynew — diag( Ml > > ■ • • , Mn ) , Mi > Mi+1 (3.20) 

This new realization {A, B, C , } is called a frequency weighted balanced realization. 
Here the balancing transformation matrix T is found using the algorithm given in [61]. 

Now the controller reduction is achieved by eliminating the rows and columns of 
A ' 13 aud C corresponding to smallest (Mr+i, Mr+ 2 , , Mn) in U new = Y new . This 

yields A,.(c). This part has been explained in Section 3.2 (see equations (3.8) - (3.9)). 
Here', unlike in nonweighted case (Section 3.2), a stable K(z) may not yield a stable 
Me). H' should be expected that concentration on stability criterion for minimization 
can lead to poorly performing controllers, while concentration on other performance 
measures for minimization could lead to instability [6]. 

3.3.1 A State Space Representation For W{z ) and V{z ) 

A state space representation for W(z) described in equation (3.14) can be found from 
block diagram shown in Figure 3.2. Let the state vectors associated with lifted plant 


Figure 3.2: Block diagram for obtaining a state space representation for W(z) 

P[z) and controller K(z) be Xi(k) and x 2 (k) respectively. State models of lifted plant 
and controller are given in Section 2.3. Let the input vector be u(k), output vector 

x i(k) 

x 2 (AO 




y(k) and state vector x(h), where x(k) 
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The state equations for lifted plant can be written as follows 


Xi(k + 1) — Ai X\(k) + Bi Ei e(k) (3.21) 

and y(k) = C l X\{k) + D l E\ e{k) (3.22) 

and the state equations for controller as 

x 2 {k+l) = Ax 2 {k) + By x {k) (3.23) 

and y 2 {k) = C x 2 (k). (3.24) 

Here 

e(k) = u(k) — y 2 (k) (3.25) 

and yi(k) = E 2 y{k). (3.26) 


Putting eqn. (3.24) in (3.25) and then putting eqns. (3.25) and (3.26) in (3.21), (3.22) 
and (3.23) respectively, we get 

'xi(* + l)l [ A l -B l E 1 C 1 r Xl (k) 1 f B L E l ] /IN 

— + u{k) 

x 2 (k + 1) J BE 2 C L A-BE 2 D l E 1 C J [ x 2 (k) J BE 2 D L E 1 _ 

1 ^l(^) 

y{k) = C L -D l E x C +[D L E x ]u(ky (3.27) 

L J [ x 2 (k) _ 

Now A w . B W ,C W . D w obtained from eqn. (3.27) are given below 

Al —BlE x C BlE x 

A w — i Bw — 

BE 2 C l a- BE 2 D l E x C J [be 2 d l e 1 _ 

C w = Cl -DlExC , D w = [D l E 1 ]. (3.28) 

Now minimal state variable realization for W (z) is found by deleting the uncontrollable 
and unobservable states from eqn. (3.27). 

A state space representation for V(z ) described in eqn. (3.14) can be found from 
block diagram shown in Figure 3.3. The state equations can be written similarly and 
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arc solved for x(k). x(k) so obtained is given below 

Ai —BiEiC xi (k) 0 

B E 2 Cl A B E 2 Dl E\ C X 2 (A;) B E 2 

... r i xi(k) 

y(k) = -E 2 C l -E 2 D l E 1 C +[E 2 ]n(k). (3.29) 

L J [ x 2 {k) _ 

Now .4,,, B r , C v , D v obtained from eqn. (3.29) are given below 

Al — Bl E\C 0 

be 2 c l a- be 2 d l e 1 c v be 2 
c v = -E 2 C l -E 2 D l E 1 C , D v = [E 2 ]. (3.30) 

Now minimal state variable realization for V(z) is found by deleting the uncontrollable 
and unobservable states from eqn. (3.29). 

3.4 Construction Of Balancing Transformation 

An algorithm [61] is presented below which computes the state space balancing trans- 
formation defined in Section 3.2 and Section 3.3. 

Step 1 : Compute Cholesky decomposition of Gramians U and Y : 

Let R c and R 0 denote the lower triangular Cholesky factors of the Gramians U 

and Y, i e- 





U = R c Y = R 0 R t 0 . 
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Step 2 : Compute singular value decomposition of the product of the Cholesky factors: 
Compute the SVD 



R X R C = G ES T 

(3.31) 


T = diag ( cp , <r 2 , . . . , a n ), op > cr 2 > ... > o n > 0. 


Step 3 : Form 

the balancing Transformation: 



T = R C S E- 1 / 2 . 



also T - 1 = E~ 1/2 G t R x . 

(3.32) 

hero 

T~ 1 UT~ t = E. 

(3.33) 

and 

T t YT = S. 

(3.34) 

From equations 

(3.31), (3.33) and (3.34) it is clear that the controllability and observ- 


ability gram inns are equal and diagonal. 


3.5 Implementation of the Controller Order Reduction Pro- 
cedure 

First, wo s oo how to obtain K for the LQG problem. The continuous time plant is 
described by equation (2.1). The discrete time model of the plant can be found from 
equations (2. 2-2. 3), for a sampling period of r secs., N = 1 and a zero order hold. The 
continuous time model of the plant is given as follows 

= Xp ^ + Bp 

z p (t) = MpX p (t) + V p u(t) + v p (t). (3.35) 

Whore w p and v p are assumed to be mutually independent, zero mean, continuous time 
noise processes with positive definite covariance matrices W c and V c , respectively. 

The continuous time controller is described by the following equation 

= A p x c (t) + B p u(t)+H c (z p (t)-M p x c (t)) 


(3.36) 
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and control input by 


u(t) = -F c x c (t). (3.37) 

Here 

H c : Kalman filter gain matrix for continuous time plant 

F t - : Feedback gain matrix which is obtained so as to minimize the following quadratic 
cost function [60] 

roo 

J c = / {zJ(t)Q c z p {t) + u T (t)Ru(t))dt (3.38) 

J 0 

where u and z v are control and measurement vectors, Q c is symmetric positive 
sc'midefinitc matrix and R is positive definite weighting matrix of appropriate 
dimensions. For the existence and uniqueness of the solution to the LQR we 
assume standard controllability and observability conditions. 

The discrete time model of the plant is given as follows 

X d (^k d~ 1) — -Trf (/c) “l - B d u(^k^j "F 'll) d 

z d (k ) = M d x d {k) + V d u{k) + v d . (3.39) 

Whore w d and v d are assumed to be mutually independent, zero mean, discrete 
noise processes with positive definite covariance matrices W d and V D , respectively. 

The corresponding discrete time controller is described by the following equation 
x(k + 1) = A d x(k) + B d u(k)+H d (z d (k)-M d x(k)) (3.40) 
and control input by 

u(k) = -F d x(k). (3.41) 


Hero 


H d : Kalman filter gain matrix for discrete time plant 
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Fd - feedback gain matrix which is obtained so as to minimize the following 
quadratic cost function [60] 

OO 

J d = E Y2i z d(k)Q d z d (k) + 2zJ(k) M m u(k) + u T (k) R d u(k)) (3.42) 

k= 1 

where u and z d are control and measurement vectors, Q d is symmetric positive 
semidefinite matrix and R d is positive definite weighting matrix of appropriate 
dimensions. For the existence and uniqueness of the solution to the LQR we 
assume standard controllability and observability conditions. J d in eqn. (3.42) is 
found by descretizing J c in eqn. (3.38) and it in general includes a cross product 
term AI m in the performance index even though J c may not contain any. The 
continuous time controller matrices A c , B c , C c are obtained as given below 

A c = A p - B p F c - H c M p 
B c = H c 

C c = -F c (3.43) 

and the discrete time controller matrices A, B, C are obtained as given below 

A — A d — B d F d — H d M d 

B = H d 

C = -F d . (3.44) 

The LQGs controller can be found by descretizing the equation (3.43) for a 
sampling period r secs, and a zero order hold. 

The state feedback gain matrices F c and F d and the Kalman filter gain matrices 
H c and H d may be determined by solving the standard LQR and LQG problems 
for continuous time and discrete time problem respectively. In this thesis these 
are obtained with the help of software package 3NLATLAB (version 3.12). 
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3.6 Conclusions 

In this Chapter, approximation of high order controller by low order controller has 
been considered. Controller order reduction by balanced realization truncation 
and frequency weighted balanced realization truncation has been presented. In 
case of frequency weighted balanced realization truncation a weighting function 
for preserving the closed loop transfer function has been presented. Also, a 
balancing transformation matrix has been presented. Performance criterion for 
designing the robust optimal controller has been presented. 



Chapter 4 


NUMERICAL EXAMPLE 


4.1 Introduction 

In this chapter, a sufficient condition for linear system stability robustness is 
presented. Also, effectiveness of the theory presented in the previous chapters is 
seen by designing the digital autopilot for the launch vehicle. 

Performance of various reduced order controllers are compared with the existing 
9 th order LQG controller (with lifting and without lifting the system). The 
element bounds in the system matrices for the allowed perturbations with the 
designed controllers are presented. The step responses of various reduced order 
controllers and 9 th order LQG controller are shown and compared with each 
other. 

4.2 A Sufficient Condition For Stability Robustness 

Consider a linear discrete-time closed loop system described by 

x(k + 1) = Ax{k). (4.1) 

It is assumed that stable nominal system (4.1) is subject to linear time-varying 
parametric uncertainties in the entries of A described by A A(k ) . The closed loop 
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system matrix is given in equation (B.l;. Thus the linear discrete-time closed 
loop systems with time- varying uncertainties is given by 

x(A--fl) = (.4-A.4(A))x(A). (4.2) 

Let bounds be given on the absolute values of the maximum variations in the 
element of A .4(A). i.e. 

A.4[} < max ( A. 4~i AT) VA = 0.1 (4.3) 

and A. 4 is defined as a matrix with entries max A. 4+(A). A. 4(A) and A .4" (A) 
matrices are given in eqns. (B.2) and (B.3) respectively. Then the structured 
uncertainties are defined as 

{ A. 4(A) : A AT (A) [ < ] A .4 V A = 0,1... }. (4.4) 

Here in A .4+ (A) [<] A«4 ‘ <' is applied eiement-by-element to two matrices. 

.4+ : A matrix obtained by replacing the entries of A with their absolute values. 
Further, we define 

e = max ( A A,., ) 1 = 1, ... ,n; (4.5) 

j = 1 - - At 


and E is a matrix with entries 

Eij = - (4-6) 

J e 

Then the structured uncertainty (4.4) can be rewritten as 

(AA(A) : AA + (k)[<]eE VA =0,1...}. (4.7) 

Now the stability robustness problem can be stated as : 

Given a structured uncertainty matrix E, determine the upper bound on e for 
maintaining asymptotic stability of the closed loop system (4.2) for ah A 1(A) 
described by (4.7). The e can be found with the help of theorem given in [59], 

see also [58]. 
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Now in the following theorem, a sufficient condition for stability robustness is 
given. 

Theorem 4.1 [59] 

Suppose the nominal closed loop system described by (4.1) is asymptotically 
stable, and M is the modal matrix to transform the nominal closed loop system 
matrix .4 to a Jorden form matrix B. Then the uncertain closed loop system 
given by (4.2) is asymptotically stable for all A A{k) described by (4.7), if 

< 7r [ (/ - B + ) _1 (M -1 )+ EM+Y 

where 

~(A) : Perron-Frobinius radius, i.e., the maximum eigenvalue of 
non-negative matrix A. 

4.3 Design Of Autopilot For Launch Vehicle 

4.3.1 System Model 

The launch vehicle has been modeled as 1-input, 3-output, 23 rd order continuous- 
time system. The discrete-time model is obtained from the corresponding continuous- 
time model for a sampling period of 0.04 secs., A = 2 and a zero order hold. Zeros 
of the full order continuous-time nominal model as well as those of the two sets 
of perturbed continuous time model are given in Table A.l. Zeros of the discrete 
time nominal model and the perturbed models are given in Table A.2. 

The reduced order model of the launch vehicle is of the 9 th order. Both the 
continuous-time and discrete-time full order as well as reduced order models are 
nonminimum. The reduced order continuous-time and discrete-time model are 
used for the design of the controller. Loop transfer recovery is only partial with 
the nonmimimum system. Controller order is further reduced using technique 
presented in Chapter 3. The designed controller is implemented on the full order 
system. 
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The reduced order continuous-time mode: is characterized by matrices in (A.l). 
(A.2j. iA.3). (A. 4). (A. 5) and (A. 6). The discrete-time model is obtained from 
the corresponding continuous-time model for a sampling period of 0.04 secs.. 
A = 1 and a zero order hold using equation (2.2). State weighting matrices Q c 
and Qd are given by (A. 7) and (A. 10) respectively. Input weighting matrices R 
and R c ; are given by (A. 9) and (A. 8) respectively Covariance matrices Ii’ c and 
\V d of the fictitious system noise process injected into the system for the loop 
recovery art 1 given by (A. 11) and (A. 12) respectively. Covariance matrices \' c and 
1 'd of the measurement noise process are given by (A. 13) and (A. 14) respectively. 
Stability factor are given in (A. 15) and (A. 16) respectively. 

The above mentioned matrices are used for designing the LOG $ (Controller ob- 
tained dftcr descretizing the continuous time controller) and LQGp (controller 
obtained for the discrete time model of the plant) controllers. The designed con- 
trollers should satisfy following performance measures for the closed loop system: 

Overshoot less than 6%. 

Gammargin > G dB. 

Phaftcnuirgtn > 35°. 

Bandwidth approximately equal to or slightly greater than 5 rad/sec.. 

Feedback should tolerate a delay of 40msec. 

The designed controllers must work satisfactorily for the two perturbed models 
of the launch vehicle. 

4.3.2 Robust Optimal Controller Design 

For simplicity following notations are used 

LQGbd Balanced realization of LQGd controller. 

LQGbs Balanced realization of LQGs controller. 
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LOG fed Frequency weighted balanced realization of LQG D controller. 

LOG fbs Frequency weighted balanced realization of LQG S controller. 

To determine .4. B and C matrices for the controller Kiz) and ,4 C . B c and C c for 
the controller K c (t). Feedback gain matrices F d and F c and the Kalman filter gain 
matrices H d and H c are determined as mentioned in Section 3.5. The feedback 
cam matrices F d and F c obtained are given in (A. 17) and (A. 19) respectively. 
The Kalman filter gain matrices H d and H c obtained are given in (A. 18) and 
A. 20) respectively. From F d and H d .4. B and C matrices for LQG D controller 
are found using equation (3.44). From F c and H c .4 C . B c and C c matrices for 
A'ctfi controller are found using equation (3.43). Now the LOG s controller is 
found by descretizing the I\ c (t) controller for a sampling period of 0.04 secs.. 
A’ = 1 and a zero order hold. 

With the above mentioned LQGe controller and LOGs controller, now bal- 
anced realization truncation procedure as described in Section 3.2 and frequency 
weighted balanced realization truncation procedure as described in Section 3.3 
t an be applied to further reduce the order of controllers. Here. LOGe and LQGs 
t out rollers are stable. For balanced realization truncation K r (z) is chosen to min- 
imize the condition given in equation (3.1). The balancing transformation matrix 
T is found by carrying out steps presented in Section 3.4. With this T matrix .4*. 
B b and C b are determined using equation (3.6). This new realization {A b . B b . C b } 
is a balanced realization. 

For frequency weighted balanced realization, weights TT*(c) and K(c) are formed 
as described in equations (3.14-3.15), and U and Y are determined from equa- 
tions (3. 17-3. IS). Now the balancing transformation matrix: T is found out by 
carrying out steps presented in Section 3.4. W ith this T matrix A b . B b and C b are 
determined using equation (3.6). This new realization {A b . B b . C 6 } is a frequency 
weighted balanced realization. 
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The A b . B 0 . C b . and T matrices obtained for LQG B d controller are given by 
iA.21). (A. 22;. (A. 23). (A.24) and (A.25 respectively. The A b . B b , C b , E and 
T matrices obtained for LQGbs controller are given by (A. 26). (A. 27). (A. 28). 
i A. -9) and .A. 30) respectively. The Aj,. B b . C b . E and T matrices obtained 
for LQG fb ? controller are given by (A.31). (A.32). (A.33). (A.34) and (A.35) 
respectively The A b . B b . C b . E and T matrices obtained for LQG FBD controller 
arc given by >. A.36). (A. 37). (A. 38). (A. 39 and .A. 40) respectively. 

Xcnv „-l r . B. and C r for balanced realization truncation and frequency weighted 
balanced realization truncation can be found as mentioned in equation (3.9). This 
constitutes the required reduced order controller. 

Performance of various reduced order controllers, with and without lifting the 
system are considered next. Lifting of plant and the closed loop system is done 
as mentioned in Section 2.2. Here r = 0.04 secs, and X = 2 is taken. So the 
plant is fast sampled at 0.02 secs, and the controller is sampled at 0.04 secs. This 
constitutes multi rate problem. Using lifting technique this multirate problem 
is replaced by a single-rate problem. The sampling frequency (314.16 rad/sec.) is 
much higher than the closed loop bandwidth (5 rad/sec.), so antialiasing filter 
A not needed in our case. Transfer function of the lifted plant is found using 
'■(juations (2.6-2. 7). 

Table 4.1 compares performance of the LQGbd • LQG B s ■ LQG FB s and LQG FBD 
controllers for different controller order. Here fast sampling of the plant is taken 
into account. As can be seen from Table 4.1 4^ order LQG B d controller, 3 rd 
order LQGbs controller, 7 th order LQG FB d controller and 2 nd order LQG FBS 
controller results in lowest % overshoot for the respective controllers. Also, % 
overshoot in case of LQG B d and LQG FB d controllers is less than in case of the 
LQGbs and LQG fbs controllers for different controller order. Also, it is seen 
that performance of 2 nd order LQG fbs controller is better than other 2 nd order 
controllers presented. This may be true for the numerical example presented 
but in general this may not be true. Also, it is seen that performance of above 
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controller 

order 

Maximum overshoot 

LQGbd 

controller 

LQGbs 

controller 

lqg fbd 

controller 

LOG bos 

controller 

9 

2.6 

7.58 

2.6 

7.6 

8 

2.6 

7.58 

2.6 

7.6 


2.6 

7.58 

2.6 

7.58 

6 

1 

7.7 

4.56 

8.87 

5 

1.8 

5.9 

4.25 

6.69 

4 

2 

9.11 

21.6 

8.3 

3 

sluggish 

4.5 

sluggish 

12.S4 

2 

oscillatory 

8.10 ! 

sluggish 

2.05 

i 

sluggish 

80 

sluggish 

13.65 


Table 4.1: 9c overshoot with lifting 

mentioned reduced order controllers are better than original 9 th order LOGp 
and LQGs controllers. This confirms the findings of other researchers that in 
many circumstances reduced order controller performs better than high order 
controller. 

Table 4.2 compares performance of LQGbd • LQGbs - LQGfbd , LQGfbs con- 
trollers for different controller order and without lifting. As can be seen from 
Table 4.1 and Table 4.2, lifting results in better performance (less overshoot) . 
Here with multirate sampling overshoots between sampling instants is reduced 
considerably compared to single rate sampling. It was seen that further lifting 
(incresing the sampling rate for the plant) did not improve the performance. 

The LQGd and LQGs controllers are designed for the reduced order model of 
order 9, but are implemented on the full order model of order 23 of the launch 
vehicle. The step responses of the state x x (pitch angle) for the 4 th order LQGbd-, 
3 rd order LQGbs-, ~ th order LQG FB d and 2 nd order LQG fbs controllers for the 
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i 

t 

. controller 

i 

order 

Maximum overshoot 

LQGbd 

controller 

LQGbs 

controller 

LQGfbd 

controller 

LOG bds 

controller 

9 

2.65 

7.65 

2.65 

7.65 

s 

2.62 

7.65 

2.64 

7.65 

7 

3.32 

7.66 

2.61 

7.65 

6 

1.3 

l . 1 

4.7 

8.95 

5 

1.8 

1 6.41 

4.4 

6.9 

| 4 

2 

9.11 

21.75 

8.4 

3 

sluggish 

4.92 

sluggish 

12.95 

2 

sluggish 

8.32 j 

sluggish 

2.15 

i i 

sluggish 

82 

sluggish 

13. S 


Table 4.2: 9c overshoot without lifting 

full order system are compared with those of the 9 th order LQGd and LQG$ 
controllers for the nominal model and two sets of perturbed models of the launch 
vehicle in Figures 4.1 — 4.6. Here step responses for 7 th order LOGfbd controller 
are shown. As can be seen from Table 4.1 controller order can be reduced to 5 th 
order if desired. 

The step response of state x.\ for 4 th order LQGbd and 7 th order LQGfbd con- 
trollers for the nominal model (Figure 4.1) shows reduction in maximum over- 
shoot compared to 9 th order LQG D controller. The step response of state ;ri for 
3 rd order LQGbs and 2 nd order LQGfbs controllers for the nominal model (Fig- 
ure 4.2) shows reduction in maximum overshoot compared to 9 th order LQG S 
controller. Comparing Figure 4.1 and Figure 4.2, it is found that rise time is less 
in case of the LQG S , LQG S s and LQGfbs controllers compared to the LQG D , 
LQGbd and LQGfbd controllers. 
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Gain margin fdB) 

Phase margin (in deg.) 

Controller 

Bandwidth 

with 

without 

with 

without 


(rad /sec.) 

lifting 

lifting 

lifting 

lifting 

9 th order LQGb 

5.05 

24.08 

23.64 

50.26 

46.56 

9 th order LOGs 

6.02 

14.06 

12.03 

49.65 

44.49 

4 tk order LQG B d 

5.12 

24.31 

23.94 

51.50 

47.71 

j 3 rd order LOG B s 

6.08 

14.09 

12.02 

49.91 

44.79 

7 th order LQGfbd 

5.07 

24.09 

23.69 

51.35 

47.51 

j 2 nd order LOGfbs 

6.12 

15.59 

14.31 

50.07 

1 

44.98 


Table 4.3: Bandwidth. Gain margin and phase margin with different controllers 

Comparing the step responses of the state Xi for 4'* order LOG BD , 3 rd order 
LOGbs- ~ th order LQGfbd and '2 nd order LQGfbs controllers for the perturbed 
system models set I and set II. with that of 9 th order LQG D and LOG s controllers 
shown in Figures 4.3 -4.6. it is observed that the peak overshoot and undershoot 
are substantially reduced in case of 4 th order LOG BD ■ 3 rd order LQG B s ■ ~ th order 
LOGfbb and 2 nd order LQGfbs controllers. From Figures 4. 3-4. 6 it is seen that 
settling time is smallest m case of the 2 nd order LOGfbs controller. From Figure 
4.5 it is seen that settling time is large in case of 4 th order LQG BB controller 
though there is reduction in maximum 91 overshoot. Taking all the factors (like, 
91 overshoot, rise time, settling time and controller order) into account it can be 
concluded that performance of 2 nd order LQGfbs controller is better than any 
other controller in two different representation presented for any order. 

Two other measures of comparison come from the gain margin and phase margin. 
Gain margin and phase margin are calculated with and without lifting the system. 
The results obtained are summerized in tabular form and is presented in Table 


4.3. 
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As can be seen from Table 4.3 gain margin and phase margin are more for 9 th order 
LQG D , 4 th order LQGbd and 7 th order LQGfbd controllers compared to 9 th 
order LQG S , 3 rd order LQG BS and 2 nd order LQG F bs controllers. Gain margin 
and phase margin are more with lifting the system for the designed controllers. 
Also there is small variation in gain margin and phase margin for high order 
LQGd and LQGs controllers and reduced order LQGq and LQGs controllers in 
balanced realization and frequency weighted balanced realization. A time delay 
of 40 msec, was put in the feedback and the closed loop system was found 
to be stable. Also bandwidth was calculated with the designed controller and 
the results obtained are given in Table 4.3. Bandwidth for 2 nd order LQGfbs 
controller is found from Bode magnitude plot shown in Figure 4.7. It can be seen 
from Figures 4. 3-4. 6 that the designed controllers work satisfactorily for the two 
perturbed models of the launch vehicle. From the above discussions it can be 
concluded that the designed controllers satisfy all the performance measures that 
have been specified for the closed loop system. 

Bounds on allowed perturbations in system matrices are determined as mentioned 
in Section 4.1. Bounds are calculated with and without lifting the plant. The 
results obtained are summerized in tabular form and is presented in Table 4.4. 

As can be seen from Table 4.4 with lifting, closed loop system is asymptotically 
stable with more perturbations in system matrices than without lifting. This 
clearly shows that with lifting, stability robustness to plant parameter perturba- 
tion is achieved. Bounds on allowed perturbations in system matrices obtained 
are compared with the actual perturbations in each entries of the system matrices 
(this is determined from two sets of perturbed models of the launch vehicle). The 
actual maximum perturbations in each entries of the closed loop system matrix 
(containing plant and the controller) for the two sets of perturbed models are as 
follows 

maximum perturbation for the perturbed system model I : 1.95 x 10~ 5 . 

maximum perturbation for the perturbed system model II : 2.60 x 10 -5 . 



38 


J Controller 

j e with lifting 

e without lifting 

9 th order LQG D 

1.019 x 10~ 4 

2.697 x 10~ 5 

9 th order LQG S 

1.140 x 10" 4 

2.741 x 10~ 5 

| 4 th order LQGbd 

1.093 x 10~ 4 

2.673 x 10“ 5 

\ border LOG B s i 1-138 x lO" 4 

2.916 x 10~ 5 

i 7 th order LOGfbd j 1.119 x 10~ 4 

3.250 x 10~ 5 

; o th order LQG F bd 

1.009 x 10~ 4 

2.984 x 10- 5 

\ 2 nd order LQGfbs 

1.102 x 10~ 4 

3.105 x 10~ 5 


Table 4.4: Perturbation bound with different controllers 

By comparison it was seen that the bounds obtained include the actual per- 
turbations in system matrices. Also it is seen from Table 4.4 that 7 th order 
LOGfbd controller accomodates larger perturbation in system matrices than 5 th 
order LOGfbd controller. This shows that the controller which gives lowest % 
overshoot can accomodate larger perturbation in system matrices. 

4.4 Conclusions 

In this chapter. A sufficient, condition for stability robustness has been presented. 
Also, robust optimal controllers have been designed for the launch vehicle. The 
design of the controllers have been done using the reduced order continuous time 
model and discrete time model of the launch vehicle. Controller order reduction 
by balanced realization truncation and frequency weighted balanced realization 
truncation method have been presented. The various reduced order controllers 
have been implemented on the full order model. 

Bounds for the allowed perturbations in system matrices and step responses of 
the various reduced order controllers have been found out and compared with 
those of the 9^ order LQG D and LQG S controllers. 



Kth Sampling instant 


Figure 4.1: Step-Responses of the 9 th order LQG D , 4 th order LQG B d and 7 th order 
LOG fed controllers on full-order(23rd) nominal system. 
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Kth Sampling Instant 

Figure 4.5: Step-Responses of the 9 th order LQG D , 4 th order LQG BD and 7 th order 
LQGfbd controllers on full-order(23rd) perturbed system set-II. 





Kth Sampling Instant 


Figure 4.6: Step-Responses of the 9 th order LQG S , 3 rd order LQGbs and 2 nd order 
LQGfbs controllers on full-order(23rd) perturbed system set-II. 








Chapter 5 


CONCLUSIONS 


In this thesis, performance of a hybrid closed loop system with and without lifting is 
considered. Here, the aim is to design a robust optimal controller. Two cases has been 
considered: 

Controller designed for discrete-time model of the plant. 

Controller obtained by descretizing the continuous-time controller. 

The controllers has been designed using the reduced order model of the launch 
vehicle and is implemented on the full order model of the launch vehicle. Controller 
order is further reduced by balanced realization truncation and frequency weighted 
balanced truncation methods. A comparison has been made of the performance of the 
system with lifting the plant and without lifting the plant. Also, a comparison has been 
made of the performance of the various reduced order controllers with the previously 
designed controller for the same system. 

The comparison clearly indicates the benefit of lifting the system. It was seen that 
increasing the sampling rate further did not improve the performance.lt was seen that 
9c overshoot for the closed loop system with LQGbd and LQG FB d controllers is less 
than with LQGbs and LQGfbs controllers for different controller order. But taking all 
the factor (like, 9c overshoot, rise time, settling time and controller order) into account 
it has been found that performance of the given hybrid closed loop system with the 
designed 2 nd order controller in frequency weighted balanced realization is better than 
any other controller in two different representations presented for any order. 
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The element bounds in the system matrices for the allowed perturbations with the 
designed controllers has been obtained. It was seen that with lifting closed loop system 
is asymptotically stable with more perturbations in system matrices than without 
lifting. 



Appendix A 

In this Appendix the data for the reduced order system, the feedback gain matri- 
ces and the Kalman filter gain matrices for continuous and discrete time systems 
along v> ith balanced controller matrices are given. 


Zeros of the 

Zeros of the 

Zeros of the 

Nominal System Model 

Perturbed Model Set I 

Perturbed Model Set II 

-1.5592 xlO 6 

-0.3022 4 0.52341 xlO 6 

0.2875 4 0. 49781 xlO 6 

0.7794 4- 1.35021 xlO 6 

-0.3022 - 0.52341 xlO 6 

0.2875 - 0.4978i xlO 6 

0.7794 - 1.35021 xlO 6 

0.6044 xlO 6 

-0.5747 xlO 6 

-41.3420 4-29.50801 

51.6977 4l9.065Si 

70.8324 

-41.3420 -29.5080i 

51.6977 -19.06581 

-56.3889 

42.5823 -4-25.64101 

-49.1475 426.9139i 

-0.4S79 448.064H 

42.5823 -25.6410i 

-49.1475 -26.91391 

-0.4S79 -48.064H 

-0.4864 4-47.9400i 

-0.4856 447.8637i 

-39.0661 

-0.4864 -47.9400i 

-0.4856 -47.86371 

29.0544 

-0.4221 438.35301 

-0.4228 438.34791 

-0.3364 430.68091 

-0.4221 -3S.3530i 

-0.4228 -38.34791 

-0.3364 -30.68091 

13.3361 

13.6440 

12.6686 

-11.5938 

-11.8496 

-10.6272 

-1.0904 4l0.SS32i 

-1.0882 410.88121 

-1.0894 4l0.SS40i 

-1.0904 -10.8S32i 

-1.0882 -10.8812i 

-1.0S94 -10.88401 

-1.0174 4 9.9535i 

-0.9952 4 9.9518i 

-0.9962 4 9.95421 

-1.0174 - 9.9535i 

-0.9952 - 9.95181 

-0.9962 - 9.9542i 


-0.5896 4 5.95861 

-0.5906 4 5.951T1 

-0.5891 - 5.96311 

-0.5896 - 5.95861 

-0.5906 - 5.95171 

-0.4531 4 4.5734i 

-0.4531 4 4.5715i 

-0.4532 4 4.5667i 

-0.4531 - 4.5734i 

-0.4531 - 4.5715i 

-0.4532 - 4.5667i 

-0.0426 

-0.0640 

-0.0617 


Table A.l: Zeros of the continuous time full order system model 


































49 


Zeros of the 

Zeros of the 

Zeros of the 

Nominal System Model 

Perturbed Model Set I 

Perturbed Model Set II 

-6.0904 

-6.0066 

-5.9702 

2.0424 -f 1.14861 

2.6081 + 1.0443i 

4.1051 

2.0424 - 1.1486i 

2.6081 - 1.04431 

-0.8S01 

-0.8542 

-0.8781 

1.7881 

0.5689 -r 0.Sl06i 

0.5702 + 0.80971 

0.5669 - 0.81201 

0.5689 - 0.81061 

0.5702 - 0.80971 

'0.5669 - 0.8120i 

0.7139 - 0.6882i 

0.7140 + 0.68811 

0.8121 - 0.57201 

0.7139 - 0.6882i 

0.7140 - 0.68811 

0.8121 - 0.57201 

1.3071 

1.3124 i 

1.2853 

0.3634 + 0.2434i 

1.1354 ! 

1.1415 -f 0.09171 

| 0.3634 - 0.2434i 

1.0764 -r 0.1693i | 

1.1415 - 0.09171. 

| 1.1058 - 0.1015i 

1.0764- 0.16931 j 

1.0471 -r 0.23691 

| 1.1058 - 0.10151 

0.9S27 + 0.2459i 1 

1.0471 - 0.23691 

! 1.0135 - 0.22091 

0.9827- 0.2459i ! 

0.9349 + 0.25981 

; 1.0135 - 0.22091 

0.9060 + 0.21761 i 

0.9349 - 0.25981 

; 0.9354 - 0.23261 

0.9060 - 0.21761 

0.8616 A 0.1760i 

1 0.9354 - 0.23261 

0.8556 -1- 0.09971 j 

0.8616 - 0.17601 

I 0.8936 — r 0.15031 

0.8556 - 0.09971 ! 

0.8044 + 0.04311 

0.8936 - 0.15031 

0.7809 

0.8044 - 0.043 li 

0.8653 

0.3213 + 0.19181 

0.4578 

0.7967 

0.3213 - 0.19181 

0.3240 

-0.1102 

-0.1179 

-0.1198 


Table A.2: Zeros of the discrete time full order system model (sampling period =0.02 
secs.) 
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System matrices are given in the following equations 


1.0655 

0.0000 

0.0004 

- 0.0001 

- 0.0004 

- 0.0005 

0.0001 

0.0002 

0.0157 

- 0.0005 

0.0003 

- 0.0005 

0.0015 

- 0.0001 

- 0.0002 

0.0000 

- 0.2842 

15.6628 

- 0.0008 

0.0001 

- 0.0001 

- 0.0004 

0.0000 

- 15.6626 

- 0.2844 

- 0.0001 

0.0010 

0.0000 

0.0003 

0.0002 

0.0009 

- 0.0020 

- 0.4107 

32.9150 

- 0.0021 

0.0006 

- 0.0002 

0.0046 

0.0003 

- 32.9152 

- 0.4112 

- 0.0009 

0.0000 

0.0000 

- 0.0001 

0.0001 

- 0.0009 

0.0006 

- 1.0064 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

- 0.0002 

0.0000 

0.0000 

0.00000 

0.0000 

0.0000 

0.0003 

0.0000 

0.0000 


- 0.0007 

- 0.0008 

0.0015 

- 0.0034 

0.0001 

0.0001 

0.0002 

- 0.0028 

- 0.0049 

0.0089 

- 0.0046 

- 0.0063 

0.0008 

— 0.0055 

- 19.7922 

20.1921 

- 20.1918 

- 19.7919 
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Input matrix 

B p = 14.5943 -1.2734 34.5310 13.7632 20.9145 17.5734 15.0813 

39.5588 40.3562 (A.2) 

Output matrix 

1.1402 1.4732 0.0962 0.1S66 -0.0854 -0.0408 -1.0777 

M v = 1.1399 0.0231 0.1062 -0.0713 -0.5909 1.1713 1.0212 

22.7444 0.4548 93.7207 200.6986 64.2873 31.1656 -22.6467 

-0.0095 -0.0380 

-0.4289 -0.6922 (A.3) 


127.7070 -306.0710 J 

r 

C p = 1.1402 1.4732 0.0962 0.1866 -0.0S54 -0.0408 -1.0777 


-0.0095 -0.0380 (A.4) 
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Input weighing matrix for continuous time plant 

R = 5. 

State weighing matrix for discrete time plant 


(A. 9) 


Qa = 


0.0543 

0.0687 

0.0015 

0.0094 

-0.0018 

-0.0036 

-0.0492 

0.06S7 

0.0869 

0.0020 

0.0119 

-0.0023 

-0.0046 

-0.0623 

0.0015 

0.0020 

0.0001 

0.0003 

-0.0001 

-0.0001 

-0.0014 

0.0094 

0.0119 

0.0003 

0.0016 

-0.0003 

-0.0006 

-0.0086 

-0.0018 

-0.0023 

-0.0001 

-0.0003 

0.0001 

0.0001 

0.0017 

-0.0036 

-0.0046 

-0.0001 

-0.0006 

0.0001 

0.0003 

0.0033 

-0.0492 

-0.0623 

-0.0014 

-0.0086 

0.0017 

0.0033 

0.0446 

0.0001 

0.0002 

- 0.0000 

0.0000 

0.0000 

- 0.0000 

-0.0001 

-0.0012 

-0.0015 

- 0.0000 

-0.0002 

0.0000 

0.0001 

0.0011 


0.0001 

-0.0012 

0.0002 

-0.0015 

- 0.0000 

- 0.0000 

0.0000 

-0.0002 

0.0000 

0.0000 

- 0.0000 

0.0001 

-0.0001 

0.0011 

0.0000 

- 0.0000 

- 0.0000 

0.0000 


(A. 10) 


Covariance matrix of the fictitious system-noise process injected into the system for 
the loop recovery for the continuous time and discrete time plant are given as follows. 

W c = 1000 * I 9>< 9 + B p * BT 


W d = B d * Bj 

Covariance matrix of the measurement noise process 


(A.11) 
(A- 12) 


V r _ 


1 * hx3 


(A. 13) 
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In = 10 * /• 


3x3 


A. 14' 


where J 3x3 is an Identity matrix of dimension (3x3). 

.1 , for the Kalman filter 


Stability factor = 


.045 , for the Feedback Control law 


Stability factor 


1.2, for the Kalman filter 

1.00175, for the Feedback Control law 


Feedback gain matrix for the discrete time model of the plant 


F d = 


0.7110 0.6094 0.0639 0.1028 0.0099 0.0150 -0.2497 


0.0000 - 0.0000 

Kalman filter gain matrix for the discrete time model of the plant 


(A. 15 i 


f A .16' 


(A. I! 


H d = 


0.0293 

0.2854 

0.0007 

0.1849 

-0.0945 

-0.0024 

-0.0002 

-0.0013 

0.0044 

0.0280 

-0.1178 

0.0031 

-0.0306 

0.0676 

0.0038 

0.0271 

0.1837 

-0.0010 

-0.0038 

0.1283 

-0.0020 

-0.0034 

0.0093 

0.0010 

0.0028 

-0.0293 

0.0002 


Feedback gain matrix for the continuous time model of the plant 


(A. 18) 




0.8470 0.7621 0.0789 -0.0107 0.0218 0.0062 -0.3095 


- 0.0000 - 0.0001 


(A- 19) 



Kalman niter gain matrix for the continuous time model of the plant 


19.1221 32.1700 13.4732 

20.8505 -15.7568 -20.5089 

-2.4111 2.5769 -9.2630 

0.6645 0.7940 44.4385 

H c = -2.0895 -40.9487 -0.5780 

-5.3863 6.0233 15.4034 

-14.5430 14.0600 -14.6495 

-1.3626 -0.8599 15.2589 

0.6439 —5.7479 —5.8580 

LQGbd controller is characterised by the following equations 


0.1674 

0.3846 

0.1731 

0.1401 

0.0807 

-0.0003 

-0.1187 

0.0299 

0.8015 

-0.1410 

-0.1618 

0.0023 

0.1660 

-0.7486 

0.3840 

0.2334 

0.1703 

-0.0038 

0.0315 

-0.0S99 

-0.1824 

0.7746 

—0.3657 

-0.0036 

-0.1724 

-0.0910 

-0.0192 

0.1053 

0.0323 

0.1913 

-0.0014 

-0.0016 

-0.0009 -0.0004 

0.0048 

-0.2364 

-0.0006 

0.0015 

0.0007 

0.0071 

0.0401 

-0.5062 

0.0004 

-0.0009 

-0.0025 

0.0025 

0.0615 

0.1336 

0.0000 

0.0004 

0.0005 

0.0005 

-0.0039 

0.0270 


(A. 20) 


- 0.0001 

0.0001 

- 0.0002 

0.0000 

0.0060 

—0.5586 

0.6155 

-0.0285 

0.0413 


-0.0001 

- 0.0000 

0.0004 

0.0000 

-0.0007 

- 0.0000 

-0.0003 

-0.0001 

0.0284 

0.0026 

0.1879 

0.0141 

0.3320 

0.0216 

0.4550 

-0.1340 

0.0644 

0.8083 


(A.21) 


o\ 



00 


0.3210 

0.2659 

- 0.0000 

0.0862 

- 0.1167 

0.0026 

- 0.288 

0.0608 

- 0.0008 

0.0449 

- 0.0521 

0.0010 

0.0025 

0.0244 

0.0025 

0.0004 

- 0.0002 

- 0.0157 

- 0.0000 

- 0.0002 

- o.ooos 

0.0008 

- 0.0010 

- 0.0001 

- 0.0001 

0.0001 

- 0.0003 


C b = I 0.4262 0.0136 -0.0075 -0.0180 -0.0112 0.0007 0.0000 
L 


0.0001 0.0000 


E = diag 


0.1894 0.0566 0.0501 0.0199 0.0070 0.0003 0.0002 


0.0000 0.0000 


(A. 22) 


(A- 23) 


(A.24) 



0.4564 

-0.2235 

0.0376 

-2.0136 

0.8830 

-0.5198 

0.8408 

0.2264 

0.1181 

-0.0379 

1.8961 

-0.7744 

0.6025 

1.3105 

-0.1139 

0.0385 

-0.7463 

-0.1171 

3.4017 

0.3046 

2.4187 

-0.0772 

0.6705 

-0.1019 

0.6624 

-1.0537 

0.6646 

-12.6232 

-0.0743 

0.0223 

0.S356 

0.1774 

2.0831 

0.1301 

0.3334 

0.2693 

-0.8618 

-0.2934 

0.4389 

1.2341 

0.5701 

— 5.6150 

0.1257 

-0.3005 

-0.1357 

-1.0109 

1.6433 

0.1392 

5.2415 

-0.0908 

0.1417 

-0.5126 

-0.3818 

2.9807 

0.1782 

3.6211 

-0.0748 

0.2645 

-0.1268 

0.0963 

1.0268 

0.4071 

-2.3098 


— 7.06 . 2 

19.6624 

17.S265 ■ 

-26.3445 

8.7742 

-8.0685 

— 44.32S4 

-0.4302 

1.8231 

-4.3888 

— 20.SS43 

-0.8836 

22.17S5 

-10.6217 

10.4689 

0.5764 

-8.7941 

-2.0090 



LQGbs controller is characterised by the following equations 



-0.2164 

-0.4147 

-0.0947 

0.0308 

-0.0164 

-0.0013 

0.0019 

0.2764 

0.2968 

-0.6454 

0.0206 

0.0115 

-0.0012 

0.0012 

0.1311 

0.3832 

0.1281 

0.2761 

-0.1850 

-0.0184 

0.0229 

0.0099 

0.0128 

0.2043 

0.6563 

0.6132 

0.0130 

-0.0141 

0.0643 

-0.1294 

0.0678 

-0.3304 

0.3727 

-0.1591 

0.1273 

-0.0028 

0.0159 

0.0148 

-0.0940 

0.0879 

0.9324 

0.1190 

-0.0026 

-0.0024 

-0.0371 

0.0271 

0.0909 

0.0045 

0.2139 

-0.0014 

0.0055 

0.0037 

-0.0217 

0.0194 

-0.0366 

0.0783 

- 0.0000 

0.0000 

- 0.0000 

- 0.0000 

0.0000 

- 0.0000 

0.0001 


- 0.0001 0.0000 
- 0.0001 0.0000 
-0.0015 0.0000 

0.0004 -0.0000 
-0.0091 0.0000 

-0.00S4 0.0001 

0.299S -0.0025 
0.2684 0.0281 

0.0003 0.0000 

0.5910 0.1105 -0.0002 

0.0836 -0.1589 -0.0003 
0.0169 0.1404 -0.0002 

0.0072 -0.0205 0.0007 

B b = -0.0002 0.0116 0.0000 

0.0010 -0.0014 -0.0004 
-0.0014 0.0033 -0.0001 

0.0003 -0.0005 -0.0001 
0.0000 - 0.0000 0.0001 


(A. 26) 


(A.27) 
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C 6 = 


0.6079 -0.0735 0.0111 0.0012 -0.0014 0.0001 -0.0000 


0.0000 - 0.0000 


(A. 28) 


I! = aiag 


0.3954 0.0872 0.0410 0.0072 0.0049 0.0011 0.0001 


0.0000 0.0000 


0.3155 

-0.7322 

0.3121 

0.4132 

— 0.6153 

-0.8412 

0.0430 

0.8230 

-0.4728 

-0.6562 

-0.1462 

—0.1575 

-0.5267 

-1.0415 

-0.5362 

-0.5047 

-0.4343 

-0.0596 


0.1483 

-0.4021 

-0.2685 

0.0447 

-0.7257 

-0.4042 

1.0707 

2.8173 

-0.9735 

1.7022 

2.4591 

-0.4199 

-0.5326 

-1.0S35 

-0.9517 

—0.9536 

0.1767 

1.6S71 


1.3179 -2.0776 
-1.5454 2.7401 

-4.2401 -0.9704 
1.9664 1.7361 

1.1095 1.3015 

0.7119 0.4974 

-1.2514 0.8551 

-2.4443 -0.3618 
-0.5358 0.7991 


(A. 29) 

18.4877 

-15.2599 

-1.4024 

1.9092 

-22.3918 

—7.5447 

10.8687 

1.9473 

-3.2847 


-16.2940 

3.0780 

10.0292 

—9.5175 

-12.1829 

—8.8567 

9.4973 

34.6733 

28.0207 

9.8945 

10.2746 

15.4748 

-20.7128 

-17.4626 

3.3894 

11.6382 

11.1771 

-1.3922 


(A. 30) 
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LQGfbs controller is characterised by the following equations 


0.0192 

-0.8255 

0.4950 

-0.0818 

-0.2195 

0.0453 

0.0274 

0.0047 

0.0259 

0.9757 

-0.1678 

-0.2874 

0.0589 

0.0346 

0.0120 

-0.2537 

0.5449 

0.5941 

-0.0908 

0.0258 

0.0089 

-0.0753 

0.1600 

—0.1535 

0.1298 

-0.6053 

0.0684 

0.0620 

0.0608 

0.1122 

0.0837 

0.1073 

0.4951 

-0.1092 

0.1095 

-0.0325 

0.0107 

-0.0043 

— 0.064S 

0.2388 

1.0117 

0.0712 

-0.0021 

-0.0029 

0.0001 

-0.0418 

0.0591 

-0.1643 

0.3609 

0.0007 

-0.0044 

-0.0118 

0.0105 

0.0276 

0.0494 

-0.0260 

-0.0003 

0.0003 

0.0016 

0.0032 

-0.0013 

-0.0217 

0.0295 


-0.0050 

0.0000 

-0.0063 

0.0000 

-0.0017 

0.0000 

-0.0102 

0.0001 

-0.0157 

0.0001 

-0.0073 

- 0.0000 

0.3179 

-0.0019 

0.0567 

0.1144 

0.0313 

0.0081 


17.1194 
9.3799 
-2.6572 
1.5383 
B b = 0.0963 

-0.0422 
-0.1116 
0.0248 
-0.0017 


7.6360 

-0.0045 

-3.2504 

-0.0082 

0.0960 

-0.0029 

0.6590 

-0.0116 

0.3932 

0.0002 

-0.0725 

-0.0026 

0.0472 

-0.0001 

-0.0687 

-0.0003 

0.0055 

0.0001 


(A. 31) 


(A.32) 
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C b = 0.0141 0.00S9 -0.0098 0.0012 0.0041 -0.0009 -0.0005 

0.0001 -0.0000 (A. 33) 

H = diag 738.7644 249.0767 91.1830 40.035 10.5598 2.4909 0.2072 

0.0240 0.0003 1 (A. 34) 

j 


0.0148 

-0.0021 

0.0133 

-0.0227 

0.0388 

-0.0183 

0.2511 

0.0014 

0.0126 

-0.0043 

0.0268 

0.0045 

0.1290 

-0.0878 

—0.0067 

-0.0432 

-0.769 

-0.0827 

-0.3081 

-0.0880 

-0.2634 

0.0004 

0.0068 

-0.0900 

-0.1446 

0.1751 

0.3003 

-0.6373 

-0.0011 

-0.0110 

0.0565 

-0.0882 

-0.0319 

-0.1437 

-0.4525 

0.0153 

-0.0553 

-0.0101 

0.0727 

70.0727 

0.0152 

—0.5387 

-0.0032 

-0.0166 

0.0449 

-0.0209 

0.0185 

0.2278 

0.3847 

-0.0147 

-0.0053 

0.0221 

-0.0079 

-0.103S 

0.0918 

-0.1145 

-0.0077 

-0.0088 

0.0064 

0.0031 

0.0267 

0.0196 

0.3431 


-1.S281 1.4531 

1.4690 -3.8822 
-0.C476 -3.3247 
-0.2966 13.4130 

2.4475 3.4566 ' (A.35) 

0.5548 5.8848 

-1.3576 -6.5323 
-0.1389 4.5064 

0.9667 -0.6756 
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LQGfbd controller is characterised by the following equations 


0.3297 

0.5776 

0.4312 

0.1129 

- 0.1462 

0.2267 

- 0.0246 

- 0.1298 

0.4007 

- 0.9380 

- 0.0851 

0.1129 

- 0.2466 

0.0324 

0.1549 

0.4766 

0.3417 

— 0.6675 

0.0594 

- 0.0357 

- 0.0157 

- 0.0172 

0.5163 

- 0.1568 

0.3885 

0.4690 

- 0.4661 

0.0168 

0.0414 

- 0.0016 

- 0.1692 

- 0.1201 

0.6802 

0.4219 

0.1229 

- 0.0678 

0.2488 

— 0.0754 

— 0.0545 

— 0.1865 

- 0.6840 

0.5309 

- 0.0499 

0.0252 

0.0767 

0.0160 

0.1643 

0.0447 

0.2009 

- 0.0004 

- 0.0010 

0.0002 

- 0.0077 

0.0148 

- 0.0477 

- 0.0972 

- 0.0002 

0.0103 

- 0.0005 

0.0062 

0.0065 

- 0.0430 

0.0359 


0.0026 - 0.0001 
- 0.0033 0.0001 

- 0.0003 0.0000 

- 0.0044 0.0001 

0.0080 - 0.0006 
0.0065 - 0.0007 
- 0.1380 0.0125 

0.4852 0.1785 

- 0.1461 0.8519 


10.6163 

9.2961 

0.0078 

- 1.4024 

0.5927 

- 0.0234 

- 2.8344 

- 2.5336 

0.0148 

2.8799 

- 1.5713 

0.0327 

- 0.9204 

0.1859 

0.0076 

- 0.1463 

0.3805 

0.0246 

0.2019 

0.3284 

0.0004 

- 0.0139 

- 0.0207 

0.0000 

- 0.0270 

0.0232 

0.0010 


( A . 36 ) 


(A.37) 



Appendix B 


In this Appendix the closed loop system matrix and the accompanying pertur- 
bation matrices are given. 


The closed loop system matrix containing plant and the controller is given as 
follows 


.4 = 


A d 


-B d F d 


(3-1) 


[H d C d A d -B d F d -H d C d \ 

The linear time varying parametric uncertainties in the entries of A is given by 
A ~A)(k) 


A A{k) = 


A .4^ —A B d F d 
H d A C d 0 


©- 2 ) 


and 


A.4 + (/c) 


AA+ -A BjF d 
H d ACj 0 


(B.3) 


63 . 



References 


Jlj P. P. Khargonekar. K. Poola and A. Tannenbaum." Robust control of linear time 
invariant plants using periodic compensation." IEEE Trans Auto. Contr., vol. AC- 
30. no. 11. pp. 1088-1096. 1985. 

12: D. G. Meyer. ’"A nev class of shift varying operators, their shift invariant equiv- 
alents. and multirate digital systems." IEEE Trans. Auto. Contr.. vol. 35. no. 4. 
pp. 42 9-433. 1990. 

[3j D. G. Meyer. "Cost translation and a lifting approach to the multirate LQG prob- 
lem." IEEE Trans Auto. Contr.. vol. 37. no. 9. pp. 1411-1415. 1992. 

[4] Astrom and Wittenmark. Computer Controlled System: theory and design, PHI 
Eagievood cliffs. N.J. 1984. 

[5] E C. Ifeachor and B. YV. Jervis. Digital Signal Processing: .4 Practical Approach. 
Addison-YYesley. vokingham. England 

[6] B. D. O. Anderson and Y. Liu, ''Controller reduction: Concepts and Ap- 
proaches." IEEE Trans. Auto. Contr., vol. 34. no. 8. pp. 802-812, 1989. 

[7] J. Sklansky and J. R. Ragazzini,” Analysis of errors in sampled data feedback 
systems-Partll,” AIEE Trans., vol. 74, pp. 65-71, 1955. 

[8] G. M. Krank, "Input Output analysis of multirate feedback systems,” IRE Trans. 
Auto. Contr., vol. AC-3, pp. 21-28, 1957. 

[9] J. R_ Ragazzini and G. F. Franklin,” Sampled Data Control System,” New 
York:McGraw Hill, 1958. 


64 . 



65 


[10_ R. E. Kalman and J. E. Bertram,” A unified approach to the theory of sampling 
systems." J. Franklin Inst., no. 267. pp. 405-436. 1959. 

[11] J- H. Davis." Stability conditions derived from spectral theory:Discrete systems 
with periodic feedback." SIAM J. Contr.. vol. AC-29, pp.880-887. 1984. 

!12[ E. I. Jury and F. J. MullinAThe analysis of sampled data control systems with 
periodically time varying sampling rates ." IRE Trans. Auto. Contr.. vol. AC- 4. 
pp. 15-21. 1959. 

[13’ R. A. Meyer and C. S. BurrusJ’A unified analysis of multirate and periodically 
time varying digital filters ." IEEE Trans. Circuits Syst. vol.CAS-22. pp. 162-168. 
1975. 

[14] B. Sz. Nagy and C. Foias,” Harmonic Analysis of operators on hilbert space." New 
York: American Elsevier. 1970. 

[15] P. P. Khargonekar and A. Tannenbaum." X on euclidian metrics and the robust 
stabilization of systems with parameter uncertainty ." IEEE Trans. Auto. Contr.. 
vol. AC-29, pp. 788-793. 1984. 

[16] A. Tannenbaum." Feedback stabilization of plants with uncertainty in the gain 
factor ." Int. J. Contr. . vol. 32. pp.1-16, 1980. 

[17] A. Tannenbaum." Invariance and system theory: Algebraic and Geometric As- 
pects." Xew 1 'ork: Spring er-ver lag, 1981. 

[18] B. Bamieh, J. B. Pearson, B. A. Francis and A. Tannenbaum,” A lifting Technique 
for linear periodic systems with applications to sampled data control,” Systems and 
control letters, vol. 17, pp. 79-88, 1991. 

[19] V. Yamam oto,” New approach to sampled data control systems: A function space 
method,” Proc. 29th CD C Honolulu, pp.1882-1887 Dec.1990. 



66 


.20: M. Araki and K. Yamamoto.’’ Multivariable multirate sampled data systems: State 
space description. Transfer characteristic, and Xyquist criterion ."IEEE Trans. 
Auto. Contr., vol. AC-31, no. 2. pp. 145-154 19S6. 

121] M. C. Berg. X. Amit and J. D. Powell "Multirate digital control system de- 
sign." IEEE Trans. Auto. Contr.. vol. 33. no. 12. pp. 1139-1150. 1988. 

122] D. P. Glasson/'A new technique for multirate digital control design and sample 
rate selection . " AIAA J. Guid. Contr., vol. 5. pp. 379-382. 1982. 

123] H. Al. Rahmani and G. F. Franklin,” A new optimal multirate control of linear pe- 
riodic and time invariant systems," IEEE Trans. Auto. Contr.. vol. 35. no. 4. pp.406- 
415. 1990. 

[24] H. Al. Rahmani and G. F. Franklin,” Linear periodic systems:Eigenvalue assign- 
ment using discrete periodic feedback," IEEE Trans. Auto. Contr.. vol. 34, pp.99- 
103. 19S9. 

125] A. B. Chammas and C. LeondesTPole placement by piecewise constant output 
feedback.” Int. J. Contr.. vol. 29. pp. 31-38, 1979. 

126] P. Colaneri. R. Scattoline and N. Schiavoni.” Stabilization of multirate sampled 
data linear systems." Automation, vol. 26. no.2. pp. 377-380, 1990. 

[27] B. A. Francis and T. T. Georgiou,” Stability theory for linear time invariant plants 
with periodic digital controllers,” IEEE Trans. Auto. Contr., vol. 33, no. 9, pp.820- 
832. 19S8. 

[28] D. Gangsaas, K. R. Bruce, J. D. Blight and Y. L. Ly,” Application of modern syn- 
thesis to aircraft controkThree case studies,” IEEE Trans. Auto. Contr., vol. AC-31, 
pp. 995-1104. 1986. 

[29] D. C. Hyland and D. S. Bernstein,” The optimal projection equations for model 
reduction and the relationships among the methods of Wilson, Skelton and 
Moor eC IEEE Trans. Auto. Contr., vol.ACX30, pp.1201-1211, 1985. 



67 


[30] A. G. Medievski and B. D. 0. Anderson.'' Sampled data Controller reduction pro- 
cedure." IEEE Trans. Auto. Contr.. vol.40. no. 11. pp. 922-926. 1995. 

[311 G. A. Latham and B. D. 0. Anderson," Frequency weighted optimal Hankel-norm 
approximation of stable transfer functions." System Contr. Lett., vol.5. pp.229, 
1985. 

;32[ S. \Y. Kim. B. D. 0. Anderson and A. G. Medievski." Error bound for trans- 
fer function order reduction using frequency weighted balanced truncation." Syst. 
Contr. Lett., vol.24. pp. 183-192. 1995. 

[33' B. C. Moore." Principal component analysis in linear systems: Controllability. 

observability and model reduction ." IEEE Trans. Auto. Contr.. vol. AC-26, no.l. 
pp. 17-32. 19S1- 

[34, L. Pernebo and L. M. Silverman." Model reduction via balanced state space rep- 
resentations." IEEE Trans. Auto. Contr.. vol. AC-27, no. 2. pp. 382-387. 1982. 

[35' S. Y. Kung and D. W. Lin." Optimal Hankel-norm model reductiomMultivariable 
systems ."IEEE Trans. Auto. Contr.. vol. AC-26, no. 4. pp.S32-852. 1981. 

[36[ K. Glover." All optimal Hankel-norm approximations of linear multivariable sys- 
tems and their L x -e rror bounds," Int. J. Contr.. vol. 39. no. 6. pp. 1115-1193. 1984. 

[37] A. Yousuff. D. A. Wagie and R. E. Skelton." Linear system approximation via 
covariance equivalent realizations." J. Math. Anal. AppL. vol. 106. no.l. pp. 91-115, 
1985. 

[38] A. Yousuff and R. E. Skelton,” Controller reduction by component cost analy- 
sis, 1 ' IEEE Trans. Auto. Contr., vol. AC-29, no. 6, pp. 520-530, 1984. 

[39] B. D. O. Anderson and R. E. Skelton,” The generation of all q-Markov cov- 
ers,” IEEE Trans. Circuits Syst., vol.CAS-35. no.4, pp.375-384, 1988. 

[40] A. Yousuff and R. E. Skelton,” A note on balanced controller reduction,” IEEE 
Trans. Auto. Contr., vol.AC-29, no.3, pp.254-257, 1984. 



68 


[41] J. A. Davis and R. E. Skelton,” Another balanced controller reduction approxima- 
tion." Syst. Contr. Lett., vol.4, no. 2. pp. 79-83, 1984. 

[42] B. D. 0. Anderson." Controller design: Moving from theory to practice." IEEE 
Contr. Syst. Magazine, pp. 16-25. Aug. 1993. 

[43] M. G. Safanov. R. Y. Chiang and D. J. N. Limebeer." Optimal Hankel model reduc- 
tion for nonminimal systems ," IEEE Trans. Auto. Contr.. vol.35. no. 4. pp.496-502. 
1990. 

[44' C. P. Therapos." Balancing transformations for unstable non-minimal linear sys- 
tems." IEEE Trans. Auto. Contr.. vol.34. no. 4. pp. 455-457. 1989. 

[45' J. Keller and B. D. 0. Anderson," A new approach to the discretization of contin- 
uous time controllers ." IEEE Trans. Auto. Contr.. vol.37. no. 2. pp. 214-223. 1992. 

[46] T. Chen and B. A. Francis." Input output stability of sampled data system." IEEE 
Traris. Auto. Contr.. vol.36. no.l, pp. 50-58, 1991. 

[47] R. K. Yedavalli." Perturbation bounds for robust stability in linear state space 
models." Int. J. Contr.. vol.42, no. 6. pp. 1507-1517. 1985. 

[4S] B. R. Barmish." Necessary and sufficient condition for quadratic stabilization of an 
uncertain linear system," Journal Of Optimization Theory and Application, vol.46. 
pp. 395-408, 1985. 

[49] I\. Zhou and P. P. Khargonekar,” Stability robustness bounds for linear state 
space models with structured uncertaint yT IEEE Trans. Auto. Contr., vol. AC-32, 
pp.621-623, 1987. 

[50] K. E. if ole, "Design of robust linear quadratic regulator,” American Control Conf., 
Pitsburgh, USA, pp.929-930, 1989. 

[51] O. I. Kosmidou,” Robust stability. and performance of system with structured and 
bounded uncertainty,” ini J. Contr., vol.52, no.2, pp.627-640, 1990. 



69 


[52_ J. Cieslik,"On possibilities of the extension of Kharitonov's stability test for inter- 
val polynomials to the discrete time case" IEEE Trans. Auto. Contr., vol. AC-32, 
no. 3. pp. 237-238. 1987. 

[53' A. Yicino/'Some results on robust stability of discrete time systems I IEEE Trans. 
Auto. Contr.. vol. 33. no. 9. pp. 844-847, 1988. 

[54] J. M. Martin and G. A. Hewer,” Smallest destabilizing perturbations for linear 
systems," Int. J. Contr.. vol. 45. no. 5, pp. 1495-1504. 1987. 

[55" K. M. Sobel. S. S. Banda and H. H. Yeh," Robust control for linear systems with 
structured state-space uncertainty." Int. J. Contr.. vol. 50, no. 5, pp. 1991-2004. 1989. 

[56 Q. L. Han and Q. B. Wu." Stability robustness analysis of discrete linear sys- 
tems." Int. J. Syst. Sci.. vol. 22, no.l, pp. 165-172. 1992. 

[57[ X. Niu, J. A. D. A. Garcia and E. Yaz,” Improved bounds for linear discrete time 
systems with structured perturbations,” IEEE Trans. Auto. Contr., vol. 37. no. 8. 
pp. 1170-1173. 1992. 

[5S[ J. S. Luo and A. Johnson,” Stability robustness of the discrete time LQG sys- 
tem under plant perturbation and noise uncertainty." Int. J. Contr., vol. 55, no. 6, 
pp. 1491-1502. 1992. 

[59] J. S. Luo. X. F. Zhang and A. Johnson,” Stability robustness in linear state space 
models with structured time varying uncertainty." Int. J. Contr., vol. 60, no. 3, 
pp. 435-450, 1994. 

[60] B. D. O. Anderson and J. B. Moore, Optimal Control: Linear Quadratic Methods, 
PHI, Englewood Clifis, New Delhi, 1991. 

[61] A. J. Laub, M. T. Heath, C. C. Paige and R. C. Ward,” Computation of system 
balancing transformations and other applications of simultaneous diagonalization 
algorithms,” IEEE Trans. Auto. Contr., vol.AC-32, no.2 pp.115-121, 1987. 





f| 122862 


Date Slip 


:122862 


This book Is to be returned on the 
date last stamped. 





